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ABSTRACT 

Uniformly rotating, homogeneous, incompressible Maclaurin spheroids that spin sufficiently rapidly 
are secularly unstable to nonaxisymmetric, bar-mode perturbations when viscosity is present. The 
intuitive explanation is that energy dissipation by viscosity can drive an unstable spheroid to a stable, 
triaxial configuration of lower energy - a Jacobi ellipsoid. But what about rapidly rotating compress- 
ible stars? Unlike incompressible stars, which contain no internal energy and therefore immediately 
liberate all the energy dissipated by viscosity, compressible stars have internal energy and can retain 
the dissipated energy as internal heat. Now compressible stars that rotate sufficiently rapidly and 
also manage to liberate this dissipated energy very quickly are known to be unstable to bar-mode 
perturbations, like their incompressible counterparts. But what is the situation for rapidly rotating 
compressible stars that have very long cooling timescales, so that all the energy dissipated by viscosity 
is retained as heat, so that the total energy of the star remains constant on a secular (viscous) evolu- 
tion timescale? Are such stars also unstable to the nonlinear growth of bar modes, or is the viscous 
heating sufficient to cause them to expand, drive down the ratio of rotational kinetic to gravitational 
potential energy T/\W\ oc R~ q , where R cq is the equatorial radius, and turn off the instability before 
it gets underway? Alternatively, if the instability still arises in such stars, at what rotation rate do 
they become unstable, and to what final state do they evolve? We provide definitive answers to these 
questions in the context of the compressible ellipsoid model for rotating stars. The results should 
serve as useful guides for numerical simulations that solve the exact Navier- Stokes equations in 3 + 1 
dimensions for rotating stars containing viscosity. 

Subject headings: Gravitation — hydrodynamics — instabilites — stars: rotation 



1. INTRODUCTION 

Rapidly spinning equilibrium stars are subject to non- 
axisymmetric rotational instabilities. An exact treat- 
ment of these instabilities exists for incompressible fluid 
configurations in Newtonian gravity (see, e.g., Chan- 
drasekhar 1969). For axisymmetric, uniformly rotating 
Maclaurin spheroids, global rotational instabilities arise 
from nonradial toroidal modes e bmv (m = ±1,±2, ...) 
when T/\W\ exceeds a certain critical value. Here cp is 
the azimuthal coordinate and T and W are the rotational 
kinetic and gravitational potential energies of the star, 
respectively. In the following discussion we will focus on 
the m = ±2 bar-mode, also known as the I = m = 2 
/-mode (/ = fundamental), which is the fastest growing 
mode when the rotation is sufficiently rapid. 

There exist two different mechanisms and correspond- 
ing timescales for bar-mode instabilities. Uniformly ro- 
tating, incompressible stars in Newtonian theory are sec- 
ularly unstable to bar- mode formation when T/\W\ = 
0.1375. This instability can grow only in the presence 
of some dissipative mechanism, like viscosity or gravi- 
tational radiation, and the growth time occurs on the 
dissipative timescale, which is usually much longer than 
the dynamical timescale of the system. By contrast, a dy- 
namical instability to bar-mode formation sets in when 
T/\W\ — 0.2738. This instability is independent of any 
dissipative mechanisms, and the growth time occurs on 
the hydrodynamical timescale of the system. Here we 
shall be interested in the secular bar-mode driven by 



viscous dissipation. 

Nonaxisymmetric secular instabilities in uniformly ro- 
tating compressible stars, e.g., polytropes, have been 
analyzed by several authors (see, e.g., Managan 1985; 
Imamura, Friedman & Durisen 1985; Ipser & Lindblom 
1990, 1991). The m = 2 bar-mode sets in again at 
T J\W\ ~ 0.14. However, this mode is reached only 
when the polytropic index of the star satisfies n < 0.808 
(James 1964). Stars with larger n are too centrally con- 
densed to support high enough spin in uniform rotation 
without undergoing mass-shedding at the equator (see 
Tassoul 1978 and Shapiro & Teukolsky 1983 for a dis- 
cussion and references). This constraint does not apply 
to differentially rotating stars, which can support signifi- 
cantly more rotational energy in equilibrium, even when 
the degree of differential rotation is only moderate. The 
critical value for the onset of the secular m = 2 bar- 
mode is again T/\W\ ~ 0.14 for a wide range of angu- 
lar momentum distributions and barotropic equations of 
state (see, e.g., Ostriker & Bodenheimer 1973; Bardeen 
et al. 1977; Friedman & Shutz 1978a, b). Hence stabil- 
ity criteria derived for uniformly rotating, incompressible 
spheroids turn out to have general applicability to rapidly 
rotating stars with more realistic density and angular ve- 
locity profiles, at least insofar as the m — 2 bar-mode is 
concerned. 

The point of onset of the secular bar-mode instabil- 
ity in homogeneous, axisymmetric Maclaurin spheroids 
exactly coincides with the point of bifurcation where the 
Maclaurin sequence branches off into the nonaxisymmct- 
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ric Jacobi sequence. Incompressible Jacobi ellipsoids are 
homogeneous, uniformly rotating, triaxial equilibria that 
have lower total energies than corresponding Maclaurin 
spheroids of the same mass, density and angular mo- 
mentum. The nonlinear, secular evolution of a slightly 
perturbed, unstable Maclaurin spheroid to a triaxial Ja- 
cobi ellipsoid was explicitly demonstrated by Press and 
Tcukolsky (1973), who integrated the full set of Riemann- 
Lebovitz ordinary differential equations for incompress- 
ible ellipsoids with viscosity. Sequences of uniformly ro- 
tating polytropes with n < 0.808 also have a point of 
bifurcation beyond which they are secularly unstable; 
for example, this point occurs at T/\W\ = 0.1298 for 
n = 0.75 (Ipser & Lindblom 1990), close to the value 
T/\W\ = 0.1375 for n = (i.e. incompressible) stars. 
Viscosity again is expected to drive unstable models 
to secularly stable, Jacobi-like ellipsoids in compressible 
stars. 

It seeems almost self-evident why spheroids beyond 
the bifurcation point are secularly unstable to bar-mode 
perturbations, since with viscosity such configurations 
can evolve to equilibrium configurations of lower energy, 
while conserving their mass and angular momentum. 
However, this intuitive explanation is not applicable to 
a rapidly rotating star in which the energy liberated by 
viscosity is not emitted by the star, but instead remains 
trapped in the form of thermal energy. In this case the 
star loses no mass, angular momentum or energy as it 
evolves via viscosity, and thus it cannot evolve to a lower- 
energy equilibrium state. This problem does not arise 
in incompressible stars, since they possess no internal 
(thermal) energy whatsoever and, hence, all the energy 
generated by viscosity is liberated instantaneously The 
rate of decrease in the total energy of an incompressible 
star is exactly the rate of viscous energy generation; this 
equality is guaranteed by the Ricmann-Lcbovitz equa- 
tions with viscosity. By contrast, compressible matter 
can possess thermal energy, and the fate of viscous energy 
dissipation in compressible stars depends on timcscalcs. 
If the cooling timescale, due, e.g., to thermal radiation, 
is shorter than the thermal heating timescale due to vis- 
cosity, then the thermal energy generated by viscosity 
will be radiated away quickly and the total energy of the 
star will decrease on the viscous dissipation timescale, 
as in the incompressible case. If, however, the cooling 
timescale is longer than the viscous heating timescale, 
there will be an increase in thermal energy, but the total 
energy of the star will be nearly constant on a viscous 
timescale. The possibility for the nonlinear growth of a 
bar-mode instability is not so obvious in the later situa- 
tion. In fact, in a compressible star in which the cooling 
timescale is long, any viscosity-generated enhancement in 
the thermal energy might cause the star to expand, in- 
creasing the equatorial radius R cq , decreasing the value 
of T/\W\ ~ l/i? eq and potentially turning off any non- 
axisymmetric mode that might be unstable initially. 

In this paper we address the issue of the secular bar- 
mode instability in rapidly rotating compressible stars 
in which the cooling timescale is either much longer, 
or much shorter, than the viscous heating timescale. 
We integrate the dynamical equations for a compress- 
ible star with viscosity using the formalism of Lai, Ra- 
sio & Shapiro (1994) (hereafter LRS 1994). This for- 
malism provides a set of ordinary differential equations 



(ODEs) for the evolution of the principal axes and other 
global parameters characterizing the star, which is ap- 
proximated as a triaxial, compressible ellipsoid with a 
polytropic equation of state and a velocity field that is 
a linear function of the coordinates. The equations are 
exact in the incompressible limit n = 0, where they de- 
scribe the exact hydrodynamical behavior of the objects 
and reduce to Ricmann-Lcbovitz equations identically 
This dynamical model for compressible ellipsoids is for- 
mally equivalent to the affine model developed by Carter 
& Luminet (1983,1985), although the formulation is quite 
different. 

We follow the evolution of rapidly rotating, compress- 
ible spheroids that lie beyond the bifurcation point and 
are given small, nonaxisymmetric perturbations at t = 0. 
We track their evolution in the two extreme opposite lim- 
its of very rapid and very slow cooling. We find that in 
both cases the configurations are secularly unstable to a 
bar-mode perturbation, but they evolve to different, sta- 
ble, compressible Jacobi ellipsoids. In the rapid-cooling 
limit, the final ellipsoid has lower energy than the initial 
spheroid, while in the slow-cooling limit, the final ellip- 
soid has the same energy, but a higher entropy than the 
spheroid. 

Our paper is organized as follows. In Section 2 we de- 
scribe the compressible ellipsoid model for rotating stars 
and extend the formalism of LRS to the case in which 
the entropy of the fluid is allowed to change with time 
due to heating and cooling. In Section 3 we integrate the 
dynamical equations with viscosity to show that, both in 
the case of very rapid and very slow cooling, a uniformly 
rotating, compressible spheroid that rotates sufficiently 
rapidly and is subjected to a small triaxial perturbation 
is unstable to the bar instability and evolves to a Jacobi 
ellipsoid. In this section we also show that, in fact, the 
final Jacobi ellipsoid, which is different in the two cases, 
can be determined a priori from the conserved quanti- 
ties without having to follow the evolutionary track. In 
Section 4 we confirm our numerical result that, in both 
cases, the onset of secular instability coincides with the 
bifurcation point. In particular, we sketch a proof of this 
result by examining the energy and free energy function- 
als of perturbed stars along the compressible Maclaurin 
sequence. Finally, in Section 5 we summarize our results 
and describe plans for future work. 

2. BASIC EQUATIONS 

The key dynamical equations for the evolution of a 
compressible fluid in the ellipsoidal approximation are 
derived in LRS (1994) and will not be repeated here. 
[See Section 2 of that paper for an ideal fluid, and Sec- 
tion 4 for a fluid with shear viscosity] They are obtained 
from the Euler-Lagrange equations derived from a La- 
grangian governing the dynamics of a compressible ellip- 
soid. The properties of compressible equilibria satisfying 
the stationary-state ellipsoidal equations are described in 
LRS (1993). The assumptions underlying the compress- 
ible ellipsoidal treatment, as well as comparisons between 
ellipsoidal model calculations and calculations employing 
the exact fluid equations, are summarized in these papers 
and the references therein. The basic assumption is that 
we assume that the surfaces of constant density can be 
represented approximately by self-similar ellipsoids. The 
geometry is then completely specified by the three prin- 
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cipal axes of the outer surface. Furthermore, we assume 
that the density profile p(m) inside each star, where m 
is the mass interior to an isodensity surface, is identical 
to that of a spherical polytrope with the same volume. 
The velocity field of the fluid in a frame at rest with 
the star's center of mass is modeled as a linear super- 
position of three components: (1) a rigid rotation of the 
ellipsoidal figure; (2) an internal fluid circulation with 
uniform vorticity; and (3) ellipsoidal expansion or con- 
traction. In our adaptation, where we restrict ourselves 
to the cases in which both the angular velocity and the 
vorticity are parallel to one of the principal axes, the 
exact hydrodynamical equations (the Euler and Navier- 
Stokes equations) are replaced by a set of ODEs for the 
time evolution of the principal axes, the angular veloc- 
ity of the ellipsoidal figure and the angular frequency of 
the internal circulation. As mentioned above, in the in- 
compressible limit (polytropic index n — 0) the solutions 
we derive represent exact solutions of the true hydro- 
dynamic equations. For n / 0, our solutions are only 
approximate, since the isodensity surfaces can no longer 
be exactly ellipsoidal, and the velocity field of the fluid 
cannot be described exactly by a linear function of the 
coordinates. We adopt the notation of LRS in the dis- 
cussion below. 

2.1. Energy Evolution Equation 

The equation of state characterizing the fluid is given 
by the polytropic law 

P = Kp r , T = 1 + l/n, (1) 

where P is the pressure, p is the density, T is the adi- 
abatic constant and n is the polytropic index. If the 
polytropic entropy parameter K is assumed to be con- 
stant in space and time, as in LRS (1994), the equation 
of state is barotropic with P = P(p). In this case, no 
internal energy equation is needed to describe the evolu- 
tion of the thermodynamic state of the star. In general, 
however, the gas may undergo heating and cooling, so 
that the parameter K will be a function of the changing 
specific entropy s, i.e., K — K(s). In this case the equa- 
tion of state is not barotropic, since P = P{p, s), and it 
is necessary to solve an energy evolution equation, which 
we will now derive. Combining eqn. with the first law 
of thermodynamics, 

Tds = d(e/p)+Pd{l/p), (2) 

where T is the temperature, and adopting the ideal gas 
relation for the internal energy density, e = P/(T — 1), 
yields an equation for K: 



Tds = 



P 



1 



-dK. 



(3) 



Henceforth we shall assume, for simplicity, that s, and 
hence K, are independent of position in the star, but we 
will allow them to vary with time. Integrating eqn. J3J 
over the entire compressible ellipsoid, yields an evolution 
equation for K: 



dlnK _ dS _ 

U 7. — — K 1 )—rr — 1 vis — A CO ol- 

at at 



(4) 



In eqn. (T) is the mass-averaged (mean) temperature 
defined by 

{T) = Jlf Tdm ' (5) 



S = sM is the total entropy, M is the total mass, and U 
is the total internal energy of the star, given by 



U 



(e/p) dm = kxKpl^M, 



(6) 



(eqn. 2.9 in LRS 1994), where ki is a dimensionless poly- 
tropic structure constant of order unity which depends 
on the polytropic index n (ki ~ n) and p c is the central 
density. As a simple example, for the case of an isentropic 
ellipsoid described by a Maxwell-Boltzmann equation of 
state, we have 



P = pk B T/p 1 K = p\ F k B T c /p, 



U = k 1 k B T c M/p J , 
S = U\nK/(T) = 



(T) = k 1 (T-l)T c , 
k " M lnK. 



(7) 
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In eqn. Q k B is Boltzmann's constant, p is the mean 
molecular weight and the subscript 'c' denotes that the 
quantity is evaluated at the center of the star. 

The quantity (T)dS/dt in eqn. is the net heating 
rate of the star, determined by the difference between 
the heating rate due to viscosity, r v ; s , and the cooling 
rate, A coo i, due to radiation losses. The viscous heating 
rate is given by 



d d x = -W, 



(8) 



where Ui are the components of the stellar velocity, o~ij 
is the viscous stress tensor and W is the viscous dissipa- 
tion rate. In eqn. (4.6) of LRS (1994) the quantity W is 
calculated analytically in terms of the principal axes of 
the ellipsoid a*, their velocities hi and the angular fre- 
quency of the internal fluid motion A (oc vorticity), and 
is proportional to the mass-averaged shear viscosity (u) 
defined by 

M = X7 / vdm - ( 9 ) 

Following Press and Teukolsky (1973) we parametrize 
(i>) according to (i>) = Cvisc^Gp) 1 / 2 ^ 2 , where R = 
(arenas) 1 / 3 is the mean radius of the ellipsoid, p = 
Mj (AttR 3 /3) is the mean density and C v i sc is a dimen- 
sionless constant. Choosing C V i SC < 1 guarantees that 
the viscous damping timescale will be longer than the 
hydrodynamical timescale of the star. 

The form of the cooling rate A coo i will depend on the 
details of the stellar microphysics, but here we shall only 
assume that A coo i = when K is held at its initial value 
K(s) = Kq = K(sq) in the absence of viscosity. Thus, 
we assume that the fluid would evolve adiabatically on 
the timescales of interest were the viscosity to be turned 
off completely. For example, if we take the initial con- 
figuration to be a cold, degenerate neutron star then we 
would have K — Kq > and no cooling or heating in the 
absence of viscosity (i.e. no emission from a zero-entropy 
gas). 

Setting K = Kq at t = 0, we shall treat the subse- 
quent energy evolution in two extreme, opposite physical 
regimes. In the "rapid-cooling" regime the cooling rate 
greatly exceeds the viscous heating rate for all K > Kq. 
In this case any thermal energy generated by viscosity 
will be radiated almost immediately, and the gas will 
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not be heated much above Kq. In the limit of arbitrar- 
ily fast cooling, eqn. |@J establishes a balance between 
heating and cooling, i.e., r v i S = A coo i, and maintains the 
value of K arbitrarily close to Kq. We thus have 



r • 

x VIS 



A 



cool ; 



K = Kq 



constant, ("rapid cooling"). 

( 10 ) 

In the "no-cooling" regime, the viscous heating rate 
greatly exceeds the cooling rate, in which case cooling 
is unimportant and eqn. Q gives 



dlxiK 
dt 



= r vis /U = -W/U>0, ("no cooling"). (11) 



In this case the energy generated by viscosity will go 
into thermal energy, which will be reflected in the secular 
increase of K above its initial value Kq. 

Combining the general energy evolution eqn with 
the dynamical equations for a compressible ellipsoid 
(LRS 1994, eqns. 4.10- 4.14) yields a conserved energy 
integral: 



Bit) 



I A coo \dt = E(Q) = constant, (12) 
Jo 



where the total energy E(t) is given by 
E{t) = U + W + T. 



(13) 



In the "rapid-cooling" regime the total energy E{t) of 
the star will decrease secularly by exactly the amount 
of energy dissipated by viscosity, according to eqns. ljTU)) 
and (|12|l . This is always the situation in incompressible 
stars, since they contain no internal energy and hence 
must lose energy whenever there is viscous dissipation. 
In the "no-cooling" regime, by contrast, the energy E(t) 
will be strictly conserved, as the differential rotational ki- 
netic energy dissipated by viscosity is entirely converted 
into internal energy and is not lost to the star. Compress- 
ible stars can evolve in either of these extreme opposite 
regimes, as well as in intermediate regimes, depending 
on the physical conditions. In all cases, the mass and 
angular momentum of a star is conserved in the presence 
of viscosity, but circulation is not. These quantities are 
all monitored in our evolution calculations. 

In reporting the results of our calculations below, it will 
be useful to define the following nondimensional ratios for 
energy and angular momentum: 



E = 



E 



GM 2 /R, 



-, J 



(GM 3 i? ) 1 /2 



(14) 



where 



, n -(l-ra)/(3-n) / M \ 



47T 



(n + l)K 



4nG 



n/(3—n) 



(15) 



is the radius of a nonrotating spherical polytrope of mass 
M, polytropic entropy parameter Kq and polytropic in- 
dex n. The dimensionless quantities 9\ and £i are the 
familiar Lane-Emden variables at the surface of a poly- 
trope (see, e.g. Chandrasekhar 1939). 

3. SECULAR EVOLUTION OF UNSTABLE SPHEROIDS 

We now determine the secular evolution of compress- 
ible ellipsoids with viscosity in the "rapid-cooling" and 




Fig. 1. — Viscous evolution of a secularly unstable, uniformly 
rotating Maclaurin spheroid with e = 0.94 to a stable, uniformly 
rotating, Jacobi ellipsoid. The configuration is a compressible poly- 
trope of index n = 1 that is given a slight nonaxisymmctric per- 
turbation at t = 0. The solid lines correspond to "rapid cooling" 
and the dotted lines to "no cooling" . Time is measured in units of 
the initial rotation period of the spheroid, P, and t he se mi-major 
axes ai are measured in units of Rq defined by eqn. 1151 . 



"no-cooling" regimes. Our initial configurations con- 
sist of compressible Maclaurin spheroids that are given 
small, but finite, nonaxisymmetric (triaxial) perturba- 
tions and evolved until they settle down to a final equi- 
librium state. Typical results for secularly unstable (but 
dynamically stable) spheroids are shown in Fig. 1. The 
parameters {at} are the semi-major axes of the configu- 
ration, with d3 along the rotation axis; in axisymmetry 
Oil =a,2. For the models plotted in this figure we took 
an n — 1 compressible Maclaurin spheroid with eccen- 
tricity e = (1 - al/al) 1 ' 2 = 0.94 and T/\W\ = 0.252. 
We increased the equatorial axis a\ slightly to induce 
a perturbation, setting a%l<i\ — 0.999 at the start of 
the evolution. The viscosity parameter was chosen to be 
C v ; s = 0.15; however, as long as the secular timescale 
is longer than the dynamical timescale (i.e. as long as 
C V is is sufficiently small), the secular (viscous) evolution 
track is independent of C v i s provided we "rescale" the 
evolution time according to t a 1/C V ; S . For both the 
"rapid-cooling" and "no-cooling" cases we find that the 
evolution proceeds quasi-statically through a sequence of 
compressible, differentially rotating, Riemann-S equilib- 
rium ellipsoids, finally arriving at a uniformly rotating, 
compressible, Jacobi ellipsoid. While the initial spheroid 
is unstable in cither cooling regime, both the nonlinear 
evolutionary track and the final Jacobi equilibrium ellip- 
soids are different in the two cases. 

The difference in evolution in the two cases can be ap- 
preciated by examining Fig. 2, which plots iscntropic 
(i.e. constant K) equilibrium curves for two sequences 
of compressible, n = 1, Maclaurin spheroids of differ- 
ent entropy, and the corresponding Jacobi ellipsoids that 
branch off from these sequences. The initial spheroid for 
the evolution plotted in Fig. 1 is indicated by the solid 
dot on the K = Kq Maclaurin curve beyond the bifurca- 
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Fig. 2. — Energy as a function of angular momentum for two 
isentropic, equilibrium sequences of compressible spheroids (solid 
lines) and the compressible Jacobi ellipsoids which branch off from 
them (dotted lines). The lower sequences have an entropy constant 
Ko, while the upper sequences have K = 1. 197-Kb- The initial 
unstable spheroid evolved in Fig 1. is indicated by the solid dot. 
For evolution with "rapid cooling" it evolves vertically downward 
to a Jacobi ellipsoid of lower energy (solid triangle) . For evolution 
with "no cooling" the star does not move in the figure, but it 
evolves to a triaxial ellipsoid on the equilibrium branch with higher 
entropy, K > Kg. Although viscosity drives the instability, the 
final states are independent of its magnitude. 



tion point. For the case with "rapid cooling", eqn. i|10[l 
applies; accordingly, M , J and K all remain constant 
with time while the total energy of the configuration de- 
creases, as the energy dissipated by viscosity escapes. 
The star thus evolves vertically downward in the figure to 
a Jacobi ellipsoid with lower energy E but with the same 
angular momentum J and entropy parameter K = Kg as 
the initial configuration (solid triangle). By contrast, for 
the case with "no cooling", eqn. (| 1 1 11 applies; thus, while 
M and J are again conserved, the entropy parameter K 
increases with time, since the energy dissipated by vis- 
cosity is retained and heats up the star. Consequently, 
the star evolves at constant J and E and does not move 
at all from its initial position in Fig. 2. However, the 
configuration does change - it evolves to the Jacobi el- 
lipsoid which bifurcates from a Maclaurin sequence with 
a higher entropy. The final Jacobi ellipsoid in this case 
has the same E and J as the initial Maclaurin spheroid, 
but has a higher entropy parameter K > Kq, For the 
example illustrated in Fig 1., K/Kq — 1.197 > 1 for 
this final ellipsoid. We note that unstable incompressible 
spheroids always evolve in the "rapid-cooling" regime, 
since they have no internal energy and cannot retain vis- 
cous dissipation as heat. 

Performing numerous simulations of this type, we con- 
clude that compressible Maclaurin spheroids are unstable 
to the secular bar-mode instability whether or not they 
reside in the "rapid-cooling" or "no-cooling" regimes. 
Moreover, we find numerically that the onset of the insta- 
bility coincides exactly with the bifurcation point along 
the compressible, isentropic Maclaurin curve, as in the 
case of incompressible spheroids. In the compressible el- 



lipsoid approximation, this point is located at e = 0.8127 
and T/\W\ = 0.138, independent of n (LRS 1993). How- 
ever, the evolutionary tracks and final Jacobi equilibrium 
configurations do depend on the cooling regime, as de- 
scribed above. 

Finally, we recall that sequences of uniformly rotating 
polytropes constructed from the exact (as opposed to the 
ellipsoidal) equations of hydrostatic equilibrium do not 
have a point of bifurcation when n > 0.808. We thus 
expect that the results found here for compressible ellip- 
soids also apply to exact, uniformly rotating polytropes, 
but only for configurations with n < 0.808. But they 
may also apply to a wide class of differentially rotating 
polytropes, provided T/\W\ > 0.14. 

3.1. Determining the Final Equilibrium State 

Knowing that an unstable, compressible Maclaurin 
spheroid evolves to a compressible ellipsoid allows us to 
determine the final equilibrium state without having to 
track the secular evolution. Once the appropriate cooling 
regime is assigned, the final state is uniquely determined 
by the location of the initial spheroid along the Maclau- 
rin equilibrium curve, as was illustrated in Fig. 2. Not 
only can we determine the size and shape of the final 
ellipsoid, and all other global parameters that character- 
ize the final equilibrium configuration, but we also can 
calculate its final thermodynamic state, i.e, its entropy 
parameter K, independent of the magnitude of the vis- 
cosity. It is the set of conserved quantities that allows 
us to determine the final state uniquely, given the initial 
model. 

First consider the determination of the final state of an 
unstable spheroid evolving in the "rapid-cooling" regime. 
The configuration evolves at fixed M, J and K = Kq. 
Thus the following dimensionless ratio is also constant: 

— (1 - n/5) ( —^—) = constant = j2£(«/(3-n^ 
Kn \M 3 RoJ 

(16) 

Here polytropic structure constant of order unity 

which depends on the polytropic index n (see LRS 1993 
eqn. 3.14 and Table 1). The term on the right-hand 
side of eqn. (|16|) is expressed in terms of the universal 
dimensionless parameters J and R each of which, like 
T /\W\, uniquely specifies a configuration along either a 
compressible Maclaurin or a Jacobi sequence, indepen- 
dent of n. These parameters were introduced by LRS 
(1993), (see their eqn (3.27)) where they are defined ac- 
cording to 

j2 =i (i -' /5 >G4)- (17) 

where R — (aici^&a) 1 / 3 . These universal parameters 
are tabulated along a Maclaurin sequence in Table 2 
and along a Jacobi sequence in Table 4 of their paper. 
The significance of the product on the right-hand side of 
eqn. ifTfijl is that (1) it uniquely specifies a model along 
either a Maclaurin or Jacobi sequence, and (2) once the 
initial Maclaurin model is specified, this product (un- 
like J or R separately) remains constant during "rapid- 
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Table 1. Endpoint Jacobi States for n = 1 



Initial Maclaurin Spheroid 



Final Jacobi Ellipsoids 





T/\W\ 


J 


E 


ci2 /ai 


0,3/0,1 


R/Ro 


E 


K/K 


0.8127* 


0.138 


0.299 


-0.382 


1.000 


0.583 


1.190 


-0.382 


1.000 










1.000 


0.583 


1.190 


-0.382 


1.000 


0.8200 


0.142 


0.307 


-0.378 


0.794 


0.515 


1.195 


-0.378 


1.000 










0.795 


0.515 


1.196 


-0.378 


1.000 


0.8400 


0.154 


0.328 


-0.366 


0.633 


0.449 


1.213 


-0.366 


1.000 










0.635 


0.450 


1.215 


-0.366 


1.004 


0.8600 


0.168 


0.353 


-0.352 


0.537 


0.402 


1.234 


-0.353 


1.000 










0.539 


0.403 


1.238 


-0.352 


1.008 


0.8800 


0.184 


0.382 


-0.336 


0.462 


0.362 


1.258 


-0.339 


1.000 










0.468 


0.366 


1.274 


-0.336 


1.029 


0.9000 


0.203 


0.418 


-0.317 


0.398 


0.324 


1.286 


-0.324 


1.000 










0.407 


0.330 


1.318 


-0.317 


1.058 


0.9200 


0.225 


0.463 


-0.295 


0.342 


0.288 


1.322 


-0.306 


1.000 










0.354 


0.296 


1.381 


-0.295 


1.107 


0.9400 


0.252 


0.524 


-0.266 


0.287 


0.250 


1.369 


-0.284 


1.000 










0.304 


0.262 


1.479 


-0.266 


1.197 


0.9529** 


0.274 


0.578 


-0.244 


0.249 


0.222 


1.410 


-0.266 


1.000 










0.276 


0.242 


1.577 


-0.244 


1.303 


0.9600 


0.288 


0.615 


-0.229 


0.232 


0.209 


1.437 


-0.256 


1.000 










0.257 


0.228 


1.649 


-0.229 


1.383 


0.9800 


0.340 


0.789 


-0.173 


0.169 


0.157 


1.561 


-0.214 


1.000 










0.206 


0.188 


2.050 


-0.173 


1.924 


0.9900 


0.381 


0.989 


-0.129 


0.128 


0.122 


1.692 


-0.180 


1.000 










0.179 


0.166 


2.647 


-0.129 


2.973 


0.9950 


0.413 


1.221 


-0.094 


0.093 


0.091 


1.844 


-0.148 


1.000 










0.157 


0.148 


3.477 


-0.094 


4.806 


0.9990 


0.458 


1.892 


-0.045 


0.050 


0.050 


2.231 


-0.095 


1.000 










0.139 


0.132 


6.956 


-0.045 


17.82 


1. 


0.5 


oo 


0. 


0. 


0. 


00 


0. 


1. 










0. 


0. 


00 


0. 


00 



a J, E, and Rq are defined in eqns. 1141 and 1151 : R, = (0,-10,20, 3) 1 / 3 

b One asterisk marks the secular instability point, two the dynamical instabilty point. 



cooling" evolution and thereby uniquely determines the 
final Jacobi model. 

Next consider how the final state of an unstable 
spheroid evolving in the "no-cooling" regime is deter- 
mined. In this case the configuration evolves at fixed M, 
J and E, while K increases. Suppose we guess the fi- 
nal value of K, setting K = Rs ucss > K . Then from 
eqn. Q15[l this guess yields a new value for R Q given by 

nj (3 — n) 



pgUCSS _ r> 



K 
K~o 



(18) 



Proceeding as in the previous case, we use i?Q UCSS in 

eqn. ^ to determine the product [J 2 -R ( ™ /(3_n)) ]guo SS - 
This parameter specifes a compressible ellipsoid. In gen- 
eral, however, the total energy of the ellipsoid will not 
be equal to the energy of the initial, unstable spheroid. 
Improved guesses of K gucss are required to guarantee en- 
ergy conservation. After a few iterations, the process 
converges to the desired higher entropy, compressible el- 
lipsoid with the same energy as the original spheroid. 

We have applied the above prescription to determine in 
the compressible ellipsoidal approximation the endpoint 
states of all unstable n = 1 spheroids beyond the point of 
bifurcation. Results are listed in Table 1. For each unsta- 
ble spheroid, the top row applies to the "rapid-cooling" 
regime, while the bottom row applies to the "no-cooling" 
regime. We note that the two endpoint states derived for 
the unstable spheroid with e = 0.94 are consistent with 
the results found by tracking the full secular evolution of 



this model (slightly perturbed), as shown in Figs 1 and 
2. To assess the dependence on the polytropic index of 
the star, we determine the endpoint states for unstable 
spheroids with a different polytropic index, n = 0.5, in 
Table 2. For evolution with "rapid cooling", stars with 
higher n (i.e. softer equations of state and greater cen- 
tral concentration) lose greater fractional energy than 
stars with the same initial eccentricity (i.e. T/|W]), but 
smaller n. For evolution with "no cooling" , stars with 
higher n experience less entropy increase (as measured 
by the increase in K) than stars with smaller n. The 
reason for the entropy dependence is that, according to 
eqn. ©, U k\K and ki ~ n, so that a higher value of 
K is necessary to heat a star to the same U for smaller 
n. Finally, we note that the expansion of an unstable 
star evolving with "no cooling" in comparison with one 
evolving with "rapid cooling" (measured by the ratio of 
their mean radii R) is larger for larger n, for the same 
initial eccentricity. 

4. POINT OF ONSET OF THE BAR-MODE INSTABILITY 

Our numerical integrations of the ellipsoidial dynam- 
ical equations with viscosity show that the point of on- 
set of the secular bar-mode instability in compressible 
spheroids coincides with the bifurcation point, where 
T/\W\ = 0.1375, independent of the polytropic index 
n and the cooling regime. LRS (1993) provided a rigor- 
ous proof of the first part of this assertion: they took the 
entropy parameter K to be strictly constant, effectively 
assuming "rapid cooling" , and proved that the point of 



Bar-Mode Instability 7 

Table 2. Endpoint Jacobi States for n = 0.5 a 



Initial Maclaurin Spheroid Final Jacobi Ellipsoids 





e b 


r/\w\ 




J 


E 


ci2 /ai 




R/R 


E 


K/K 





.8127* 


0.138 





.299 


-0.450 


1.000 


0.583 


1.072 


-0.450 


1.000 














1.000 


0.583 


1.072 


-0.450 


1.000 





.8200 


0.142 





.306 


-0.446 


0.796 


0.516 


1.074 


-0.446 


1.000 














0.796 


0.516 


1.074 


-0.446 


1.000 





.8400 


0.154 


0. 


.325 


-0.435 


0.637 


0.451 


1.080 


-0.436 


1.000 














0.638 


0.451 


1.082 


-0.435 


1.007 





.8600 


0.168 


0. 


.347 


-0.423 


0.541 


0.405 


1.087 


-0.423 


1.000 














0.542 


0.405 


1.088 


-0.423 


1.003 





.8800 


0.184 


0. 


.373 


-0.408 


0.469 


0.366 


1.095 


-0.412 


1.000 














0.474 


0.369 


1.107 


-0.408 


1.057 





.9000 


0.203 





.403 


-0.391 


0.406 


0.329 


1.105 


-0.397 


1.000 














0.415 


0.334 


1.128 


-0.391 


1.115 





.9200 


0.225 


0. 


.441 


-0.370 


0.350 


0.293 


1.116 


-0.381 


1.000 














0.362 


0.302 


1.158 


-0.370 


1.217 





.9400 


0.252 


0. 


.489 


-0.344 


0.297 


0.257 


1.130 


-0.361 


1.000 














0.315 


0.269 


1.206 


-0.344 


1.411 





.9529** 


0.274 


0. 


.531 


-0.322 


0.264 


0.233 


1.142 


-0.346 


1.000 














0.285 


0.248 


1.254 


-0.322 


1.647 





.9600 


0.288 


0. 


.559 


-0.308 


0.244 


0.218 


1.150 


-0.336 


1.000 














0.268 


0.236 


1.290 


-0.308 


1.849 





.9800 


0.340 


0. 


.680 


-0.253 


0.185 


0.171 


1.183 


-0.299 


1.000 














0.220 


0.200 


1.482 


-0.253 


3.377 





.9900 


0.381 





.805 


-0.206 


0.144 


0.137 


1.215 


-0.267 


1.000 














0.192 


0.176 


1.742 


-0.206 


7.058 





.9950 


0.413 





.935 


-0.167 


0.118 


0.113 


1.246 


-0.241 


1.000 














0.176 


0.163 


2.083 


-0.167 


16.45 





.9990 


0.458 


1. 


.254 


-0.104 


0.074 


0.073 


1.316 


-0.190 


1.000 














0.151 


0.142 


3.188 


-0.104 


128.6 


1. 




0.5 




oo 


0. 


0. 


0. 


oo 


0. 


1. 














0. 


0. 


oo 


0. 


oo 



a J, E, and Ro are defined in eqns. 1141 and 1151 : R, = (ait^as) 1 / 3 

b One asterisk marks the secular instability point, two the dynamical instabilty point. 



secular instability coincides with the bifurcation point for 
all n (see Section 6.2 of LRS 1993). We now sketch how 
this proof can be generalized to include the "no-cooling" 
regime as well. 

The proof of LRS is based on an energy variational 
principle and an energy functional for compressible ellip- 
soids. The functional is of the form 

E = E(p c ,X 1 ,X 2 ;M, J,K) = U + W + T, (19) 

where U , W and T are the internal, gravitational po- 
tential and rotational kinetic energies of a compressible 
ellipsoid, not necessarily in equilibrium. The ellipsoid 
of index n is parametrized by the central density p c 
and the two oblateness parameters Ai = (03 /ai) 2 / 3 and 
A 2 = (03 /CI2) 2 / 3 . The equilibrium configuration of fixed 
M, J and K (i.e. fixed entropy) is determined by ex- 
tremizing the energy according to 
f)E 

— = 0, i = 1 - 3, (20) 
don 

where {ai} = {p c , Ai, A2}. Solving these equations si- 
multaneously yields equilibrium relations for compress- 
ible Jacobi ellipsoids. Stability requires that an equilib- 
rium configuration correspond to a true minimum of the 
total energy, that is, that all eigenvalues of the matrix 
(<9 2 'E I dondoij) eq be positive. The onset of instability 
along any one-parameter sequence of equilibrium config- 
urations can be determined from the condition 

dct G3y cq =°' "- 1 - 3 - < 2i) 



When this condition is first satisfied along the sequence, 
one of the eigenvalues must change sign. By using the en- 
ergy functional for a compressible Jacobi ellipsoid given 
by eqn. I|19l) . but evaluating the determinant along the 
equilibrium Maclaurin sequence (for which Ai = A2), 
LRS found that compressible spheroids become secularly 
unstable to triaxial deformations at the bifurcation point, 
where T J\W\ = 0.1375, independent of n. 

Now suppose we allow for variations in which the en- 
tropy and entropy parameter K is not fixed, but is al- 
lowed to increase, as is the case for evolution with "no 
cooling" . In this situation, the relevant functional that 
must be minimized in a variational principle is the effec- 
tive free energy of the configuration, 



F = E — (T)S, 



(22) 



where the entropy, or equivalently, the entropy parame- 
ter K, must now be treated as another parameter that is 
allowed to vary: E — E(p c , Ai, \2,K; AI, J). The param- 
eter K appears explicitly in E only through the internal 
energy U, given by eqn. 10. However, eqn. Q shows 
that the free energy F does not change when K changes, 
while holding the other parameters fixed: 



dF_ = d^_ 0S_ = U_ dlnX = 
dK dK [ 'dK K dK 



Similarly, 



d 2 F 

doijdK 



d 2 F 
dK 2 



0, 



(23) 



(24) 
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Hence the set of parameters {o^} found to extremize E 
and identify its local minimum when holding K fixed is 
the same set which extremizes F and identifies its local 
minimum as K is allowed to vary. This set of parameters 
determines the onset of instability. Thus, compressible 
spheroids become secularly unstable to triaxial deforma- 
tions at the bifurcation point, where T/\W\ — 0.1375, 
independent of n and independent of the degree of cool- 
ing. 

5. SUMMARY AND FUTURE WORK 

We have analyzed the secular bar instability driven 
by viscosity in rapidly rotating, Newtonian stars. We 
have adopted the formalism of LRS, who model rotat- 
ing stars as compressible ellipsoids governed by a poly- 
tropic equation of state. However, we have extended 
their treatment to incorporate changes in the fluid en- 
tropy as the star evolves: the gas is allowed to heat up 
by viscous dissipation and cool down by emission. We 
treated the stars in two extreme opposite limits: one in 
which the cooling timescale is very much shorter than 
the heating timescale ("rapid cooling") and the other 
in which the timescale inequality is reversed ("no cool- 
ing"). Our numerical integrations of the ellipsoidal evo- 
lution equations show that, in both limits, a uniformly 
rotating star spinning sufficiently rapidly is secularly un- 
stable to nonaxisymmetric perturbations to a bar. The 
point of onset of the bar instability along an equilibrium 
sequence of uniformly rotating, compressible Maclaurin 
spheroids, parametrized by T /\W\, occurs precisely at 
the bifurcation point (T/|W| = 0.1375) in both regimes. 
This numerical finding is confirmed by a stability anal- 
ysis based on minimizing the free energy functional of 
the star. However, the final equilibrium state of an un- 
stable spheroid is different in the two regimes: for the 
"rapid-cooling" regime, the final configuration is a com- 
pressible Jacobi ellipsoid with the a lower energy but the 
same entropy as the initial unstable spheroid, while for 
the "no-cooling" regime, the final state is a Jacobi ellip- 
soid with the same energy but a higher entropy than the 
spheroid. The final states can be calculated from the con- 
served quantities without having to track the evolution 
of the unstable stars. 

As described in the Introduction and in LRS 
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(1993,1994), rotating compressible ellipsoids mimic many 
of the same equilibrium and stability properties of ro- 
tating stellar models that obey the exact hydrodynamic 
equations. Not surprisingly, this concordance is particu- 
larly close for stars that are governed by stiff equations of 
state (e.g. small polytropic indicies n) and are not very 
centrally condensed, since ellipsoids provide exact solu- 
tions for homogeneous, incompressible Newtonian stars. 
It is therefore relevant to examine the degree to which the 
new results found in this analysis and summarized above 
also apply to exact compressible stellar models. We are 
preparing to study the secular bar-mode instability in 
uniformly and differentially rotating, compressible stars 
by performing numerical simulations that solve the ex- 
act Navier-Stokes equations with viscosity. We will treat 
the problem both in Newtonian gravity and in full gen- 
eral relativity. We will follow not only the onset but also 
the nonlinear growth of the instability and evolve unsta- 
ble configurations until they arrive at a final equilibrium 
state. A preliminary numerical study that illustrates the 
numerical approach of tracking secular (viscous) evolu- 
tion by performing full hydrodynamical simulations is 
given in Duez et al. (2004), where we followed the secular 
evolution and determined the final fate of rapidly differ- 
entially rotating, "hypermassive" neutron stars in gen- 
eral relativity, both in axisymmetry and in 3 + 1 dimen- 
sions. Employing a hydrodynamical code to follow sec- 
ular evolution is necessary, either when the quasi-static 
configurations traversed along the secular evolutionary 
track are difficult to construct by solving the hydrostatic 
equilibrium equations (e.g., due to complicated internal 
flow patterns), or when the secular evolution ultimately 
sends the fluid into rapid hydrodynamical motion (e.g., 
catastrophic collapse). Our experience suggests that the 
technique, though resource intensive, is quite promising, 
even in general relativity. 
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